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ABSTRACT 

The dense molecular clump at the center of the Barnard 59 (B59) complex is 
the only region in the Pipe Nebula that has formed a small, stellar cluster. The 
previous analysis of a high resolution near-IR dust extinction map revealed that 
the nuclear region in B59 is a massive, mostly quiescent clump of 18.9 M . The 
clump shows a monolithic profile, possibly indicating that the clump is on the 
way to collapse, with no evident fragmentation that could lead to another group 
of star systems. In this paper we present new analysis that compares the dust 
extinction map with a new dust emission radio-continuum map of higher spatial 
resolution. We confirm that the clump does not show any significant evidence 
for prestellar fragmentation at scales smaller than those probed previously. 

Subject headings: stars: formation — ISM: clouds — radio continuum: ISM - 
infrared: ISM — radio lines: ISM 
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Introduction 



The formation of a stellar cluster proceeds when a molecular cloud clump hosts multi- 
ple dense fragments capable of collapsing independently. Cluster-for ming clumps populate 



the high-end bins of the mass and density distributions in a cloud (jWilliams et al.l 11995 



Di Francesco et al.l l2010h . but very few details are known about the process of fragmenta- 
tion, since cluster-forming clumps in very early stages of evolution are relatively rare. 



The recently popular Pipe Nebula (see lAlves et al.l 12008k and references therein), has 
become a prototype case for a cloud in a very early stage of evolution. The Pipe Nebula hosts 
only one star cluster-formin g clump, Barnard 59 (B59), and one core hosting a single young 
source in a nearby filament (IForbrich et al.ll2009f ). The rest of the cloud contain s more than 



130 starless cores that resem ble stars in the way their masses are dis tributed (lAlves et al. 



2007 ; Rathborne et al. 2009 ) and how they are distributed spatially (jRoman-Zuniga et al. 



20101 ). Dense cores in the Pipe appear to be still mostly quiescent and stable against col- 



lapse (ILada et a 



2010; Frauetal 



2008h . despite some of them being already chemically evolved ( Frau et al. 
201 la. b). Currently, the best candidate for a mechanism that supports 



(Alves et al. 


2008; 


Franco et al. 


2010) 



The B59 complex hosts one of the less massive and less distant (d = 130 pc) young 
stellar clusters observable. During the last 2.6 Myr, B59 has formed 14 stars, all below 



3 M (jCovey et al.ll2010l . hereafter CLR10). The analysis of a high-reso lution (24") near- 
infrared dust extinction map of the B59 region (jRoman-Zuniga et al. 2009 . hereafter RLA09) 
revealed that B59 is a complex group of dense cores and filamentary structures, in which 
the central clump, B59-09, hosts most of the cluster members. The analysis of the dust 
extinction map suggests that the central clump has a smooth profile that compares well with 
that of an isothermal sphere, with no evidence of internal fragmentation. Moreover, pointed 
NH 3 observations inside the core show that the thermal to non-thermal kinetic energy ratio 
averages well over unity, suggesting that B59-09 remains mostly quiescent despite having 
formed a star cluster. The main goal of this study is to confirm such a scenario, by showing 
that a radio continuum dust emission map with higher spatial resolution, also lacks evidence 
for fragmentation in the central core. 

The paper is divided as follows: in §2] we describe the observations and data reduction 
process. Our results are detailed in §3] and, to conclude, we elaborate a discussion about 
their significance in §H] . 
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2. Observations 
2.1. MAMBO-II 

We mapped B59 at 1.2 mm (250 GHz) with the MAMBO-II bolometer at the 30m IRAM 
telescope atop Sierra Nevada (Spain). MAMBO-II features a 117-receiver array that covers 
240" in diameter. The observations were carried out in November 2009 in the framework of a 
flexible observing pool. The weather conditions were good and zenith optical depth remained 
within 0.1 to 0.3. A total of 5 usable on-the-fly maps were completed and combined. The 
beam size of the telescope is 11" at the effective frequency of 250 GHz; we used a constant 
scanning at a speed of 8" s -1 in the azimuthal direction for up to 65s. This strategy resulted 
in average integration times per map of ~1 hour. Each map was performed with a different 
secondary chopping, which varied between 30" and 80", parallel to the scanning direction 
of the telescope. This procedure assured that we had a different OFF position for each ON 
position within the map. The scanning direction was also changed for each map, giving 
a different zero emission level, which helped to avoid spatial filtering effects. The zenith 
optical depth was measured with a skydip at the start and end of a map. Pointing and focus 
were checked also at the start and end of a map, with corrections below 3" and 0.2 mm, 
respectively. Flux density calibrators were observed every few hours. The data reduction was 
done with standardized routines from MOPSIC software included in the GILDAS packag^]. 



2.2. Near Infrared Dust Extinction Map and additional radio data 

We make use of the dust e xtinction map of RLA09 constructed with the NICER tech- 



nique (iLombardi &: Alvesll200l[ ) and a combined photometric catalog obtained from ground 
(ESO) and space based (Spitzer) observations. In addition to the infrared data, we also 
make use of a series of pointed observations made with the Green Bank 100m telescope to 
determine the variation of the emission of the (1,1) and (2,2) rotation-inversion transitions of 
ammonia (NH 3 ) across the central clump B59-09. These data were also used and described 
in RLA09. Finally, we also make use of a C 18 (2-1) line emission map obtained with the 
detector HERA at the IRAM 30m telescope (Roman-Zuniga et a!., in prep). 



HILDAS and MOPSIC are available at 



http://www.iram.fr/IRAMFR/GILDAS 
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3. Data Analysis and Results 

In the top panel of Figure [TJ we show our MAMBO-II map toward B59. The map 
detected mostly the emission from the clump 09at@ and, less prominently, the core 09c, at 
the northwestern end. Five young stellar objects were detected with very h igh signal-to- 



noise ratios (SNRs). They correspond to sources 6, 7, 9, 10 and 11 in the list of iBrooke et al. 



( 120071 ) (BHB07). The 250 GHz continuum emission properties of these sources are listed in 
Table [TJ Notice that sources BHB07-6 and BHB07-7 lie very close together and they are 
not resolved as separated sources in the MAMBO-II map. As our main purpose is to study 
the emission of the core, we subtracted out the contribution of the young stars. For this 
purpose, a 2D Gaussian profile was fit to each of the sources, using a background emission 
level corresponding to the average in the vicinity region of the clump. The fit parameters are 
listed in Table [TJ The smoothness of the resulting map (see central panel in Figure [lj seems 
to confirm that most of the subtracted emission arises from the YSO warm circumstellar 
material rather than from the core cold dust. 

After subtracting the contribution from the YSOs, we transformed the dust emission 
maps to column density and then into visual extinction for a direct comparison with the 
NICER map of RLA09 . The conv e rsion was done following the method described in section 



4.1 and appendix A of iFrau et al.l ( 120101 ) . We notice that the co nversion is very s ensitive to 



the assumed value for gas temperature, T k (see Appendix 1 of Frau et al. 201ol), and also 



the du st emissivity (k u ) and the dust temperature assumed. Following lOssenkopf fc Henning 



(119941 ). we assumed k„=0.007 cm 2 g 1 for grains in a dense medium (n~ 10 5 cm 3 ). Then, we 



computed visual extinction maps for gas/dust temperat ures in the 10-12 K ran ge in steps of 



0.25 K (ammonia observations yield 7jf=11.3 K±0.7 K; lRathborne et al.ll2009f ). We selected 
a temperature of 10.25 K as this value showed the best agreement with respect to the previous 
extinction map, and we made the assumption that this temperature is constant within the 
entire clump. The derived radius, mass and average density of the clump are listed in Table 
|2j Given the high SNR achieved, we restricted analysis to the region satisfying I u >0.2 I v>max 
(~4.5<r). The total mass of the clump 09ab in B59 estimated from the MAMBO-II map 
is about half of that estimated from the NICER map (i.e., cores 09ab and 09c sum about 
21 M Q according to RLA09). This difference is mostly because MAMBO-II maps are not 
as sensitive as the near-infrared excess method in detecting the diffuse gas at the external 
parts of dense cores (Ay<25 mag in this case), rather than a technical effect. The partial 
detection of cores 09cd and 04a seem to confirm the good quality of the emission maps 
despite their lack of sensitivity. The NICER technique relies on having enough sources per 



2 we follow the nomenclature of RLA09, see their Table 1 
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beam to average reddening, and background sources toward the edges of B59 are abundant. 
On the other hand, the continuum emission relies on the amount of dust that contributes 



to the signal, which decreases at lower column densities (see lFrau et al.ll2010l . for additional 
discussions on this effect). Also, the assumption of a constant temperature for the entire 
region is less accurate toward the core boundaries that might be heated by external sources. 

Morphologically, the MAMBO-II dust extinction map shows features that are equivalent 
features to those found in the NICER map of RLA09, which is shown in the bottom panel of 
FigdJ Both maps show a relatively flat central region toward the clump 09ab with a shallow 
"dent" near the location of source BHB07-10, discussed in section 5.1 of RLA09. A radial 
profile of the MAMBO-II map was constructed by averaging flux in circular, concentric rings 
centered on J2000 {a, 5) = (17 : 11 : 23,0, -27 : 25 : 59.3), as in RLA09. Figure 1 shows 
that observed profiles compare well to each other and trace equivalent structures within 
~10 AU (Ay>35 mag). Below Ay pa 35 mag, the MAMBO profile plunges and appears to 
be truncated below an average radius of 2xl0 4 AU, while the NICER profile continues until 
it reaches the background level near 5xl0 4 AU. 

We used our pointed NH 3 observations (RLA09) and the C 18 0(2-l) map to categorize 
the variation of the LSR velocity near the center of the clump. Twelve NH3 pointed obser- 
vations were made within the central 10 4 AU of B59. All of them indicate variations of vlsr 
smaller than 0.088 km s _1 from the central value of 3.485 km s _1 . Moreover, these variations 
are below the average velocity dispersion (cr v ) = 0.210 ± 0.0498 km s" 1 and below the sound 
speed in a 10 K gas (c s = 0.12 km s _1 ). The C 18 (2-1) map reveals variations below 0.05 
km s _1 from the central velocity value. The linewidths, although being significantly wider 
than those of the NH3, also show very small variations (less than 0.08 km s" 1 ) from the 
averag e value of 0.87 km s -1 , indicating very uniform kinematics near the center of B59- 



09ab. iRathborne et al.l (120091 ). discussed how kinematically independent cores in the Pipe 
Nebula must show radial velocity differences larger than c s . Following this prescription, the 
differences we observe are too small to suggest the presence of kinematically independent 
substructure. 



4. Discussion 

The radio continuum continuum map has twice the resolution of the NICER map (11" 
FWHM beam vs. 24" FWHM Gaussian filter), resolving a projected size of 1420 AU (as- 
suming a distance of 130 pc to the Pipe Nebula). We could safely say that the MAMBO-II 
map should resolve structures with minimum sizes of 2000 AU to 5000 AU. Our map, how- 
ever, does not reveal any significant substructure (other than the sources or the dent) below 
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1.5 x 10 AU. This result is in good agreement with those of iRoman-Zuniga et al.l ( 120101 ). 



who reported a possible limiting scale of fragmentation in th e Pipe Nebula of abo ut 1.4 x 10 



4 



AU. As also noted there, this scale agrees with the results of ISchnee et al.l (120 lOf ) who found 
no evidence of fragmentation in cores of the Perseus Molecular Cloud at scales of 10 3 — 10 4 
AU. 

The column density estimations made from near-IR extinction and dust emission are 
very similar near the center. Both maps fail to show evidence of substructure suggestive of 
significant fragmentation in R09ab. The dust emission map confirms that the dent is real 
and thus it suggests that source BHB07-10 is possibly affecting the material surrounding 
it. The decrement of column density near source BHB07-10, however, is equivalent to a 
decrement of only 3 percent in the total mass of the clump (RLA09). Also, CLR10 showed 
that BHB07-10 is not embedded at the level at which it is projected against the map. In 
contrast, source BH B07-11 has been suggested to be the origin of a moderate gas outflow 



(lOnishi et al.lll999l ) that seems to be carving the northern part of the clump, forming core 
09c. Thus, the feedback effect near BHB07-10 may not account for a bonafide fragmentation 
process as in core 09c. 

The young stellar cluster in B59 is likely too small to hold itself together dynami- 



cally. Following lAdams &: Myersl (120011 ) . for the B59 cluster to remain stable, its relaxation 
timescale, t r i x , should be at least larger than its formation timescale, which in turn should 
be comparable to the age of the cluster. The age of the B59 cluster is estimated as being 2-3 
Myr (CLR10). The relaxation parameter, Q r i x , i.e., the number of crossings per relaxation 
time, can be estimated in terms of the star forming efficiency, e, and the number of stars, iV*, 
approximately as N±/(10e 2 In (iV*/e)). Using the present day value, e = 0.3 (CLR10) and, 
for 10 sources projected within B59-09ab, we obtain Q rix = 3.2. Thus, given the crossing 
time, of the clump t cr rj 0.06 Myr, and t r \ x = Q r i x t cr ~ 0.2 Myr, which is much shorter than 
its formation time. It is thus very unlikely that the B59 stellar cluster will remain together 
longer than a few more crossing times. 

Exper imental fits of isothermal sphere profiles, particularly Bonnor-Ebert and Dapp- 



Basu 120091 to the NICER profile suggest that the core is already out of equilibrium, although 
it is not possible to estimate its state of evolution toward collapse (see also RLA09). The 
velocity dispersion from the NH3 data and the parameters from Table [2] yield a virial pa- 
rameter value a V i r = 0.25, which suggests that the B59-09ab clump is gravitationally bound. 
Thus, we should expect the clump to be presently moving towards collapse. As shown in 
RLA09, B59-09ab remains mostly quiescent. Moreover, our data show that the clump has 
been able to amass between 9 M Q and 19 M of gas (the latter if we consider the whole dust 
extinction structure) in an apparently monolithic structure, and we know from the age of 
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the cluster that the period of mass aggregation could be as long as about 6 crossing times 
long (CRL10). What is the mechanism that held the clump together, retarding collapse and 
fragmentation for a relatively long period? 

Feedback in the form of outflows could serve as a turbulent energy injection mechanism 
that could provide the non-thermal support required to maintain the clump against collapse. 
Several authors have provided evi dence of one and possibly two outflows from embedde d 
sources BHB07-09 and BHB07-11 JOnishi et al.lll999l : iBrooke et al.l [2OO7I : iRiaz et al.ll200gh . 
and RLA09 discussed how these could be carving structures at the outer regions of the 
clump. The presence of at least 4 YSOs located at projected distances larger than the 
Jeans length of the core (see below) could suggest that B59-09ab is more likely a remnant 
of dense gas after an episode of multiple source formation within a much larger structure. 
Even source BHB07-10 could a lso be affecting the core at a shallow level, as disc ussed above. 
The numerical experiments of iKrumholz et al.l (120071 ) and lOfiher et al.l (120101 ) have shown 
that radiative feedback can be a strong agent against fragmentation. For instance, feedback 
from protostars can inhibit the formation of binaries in a fragmenting disc scenario. These 
models, however, require stars with masses above 3 M and the mechanisms work best at 
spatial scales about one order of magnitude smaller than the observed size of B59. Since 
sources in B59 have too low masses, and are not numerous, it is thus very unlikely that they 
can either or provide non-thermal support against collapse neith er provide enough feedback 
energy to inhibit fragmentation (see, e.g.. lLongmore et al.ll2010l ). 



The other plausible candidate for a mechanism that can prevent further collapse and 
fragmentation of cores in B59 (as we ll as other region s of th e Pipe Nebula) is the magnetic 
field. Numerical studies like those of iNakamura fc Lil (120111 ) concluded that for an initially 
magnetically-subcritical cloud, a strong magnetic field is able to slow down gra yitational col- 
l apse and fragmentation, decreasing the star formation rate significantly. Also, iPrice fc Bate 
(120071 ) showed that support by magnetic fields may deter density perturbations and the 
fragmentation of discs in the case o f binar y formation. The optical polarization studies of 
Alves et al.l (120081 ) and Franco et a l. (2010) strongly suggest that the Pipe is permeated by a 
magnetic field. Moreover, iFrau et al.l (120101 ) have shown that chemically evolved cores in the 
Pipe are possibly associated with a strong magnetic field, which is suggestive of significant 
magnetic support. 



Following the formulation of iMouschoviasI (119911 ). we find a Jeans length scale of Xr,cr = 
0.12 pc, which is about twice as large as the radius of the clump in the dust emission map 
within which have demonstrated the monolithic behavior of the clump (althou gh this length 



i s very close to the radius estimated from the dust extinction map). The data of iFranco et al. 



(120101 ) indicates that B59 may be sub-Alfvenic, so in the overall region the magnetic field is 
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dynamically more important than the turbulence. 



From the Jeans length scale a nd the critical magnetic length scale, Xm = 0.91(5/ [/iG])(lx 



10 3 [cm" ]/n) ( Mouschoviasl 1 1 99 ll ) . and considering both the visual extinction and dust emis- 
sion map derived values of the clump radius, we can infer that a magnetic field strength in 
the range of 0.1-0.2 mG would be enough to support the clump. These values are larger than 
those calculated from optical polarization data of the diffuse surrounding gas. However, mag- 
netic fi^ld_js_exp_e^ted_Jo_jtren^Hi toward denser regions as the collapse process evolves 
(e.g. iFiedler fc Mouschoviaj_ll993h a nd these estimates ar e not unreasonable compared to 



other dense cores (see 



Crutcher 



19991 : ICrutcher et al.ll2004f ). 



While the profile of B59 suggests the core is out of equilibrium (Roman-Ziiniga et al., in 
prep), our continuum map shows that it has not finished collapse for a time comparable to the 
age of the stellar cluster. Magnetic field support (or an equivalent combination of supporting 
mechanisms) could have been active during such a timescale. The age of the cluster suggests 
that B59-09ab has survived for a period longer than 10 tff without completely collapsing 
or fragmenting, and it is unlikely tha t it w ill fail to evolve toward protostellar collapse. 
Experiments by iGalvan-Madrid et al.l (120071 ) suggest that pre-stellar core survival can be 
assured for at most 3-10 tff for cores with densities above 10 5 cm -3 , independently of mass 
to magnetic flux ratio. Despite the lack of fragmentation in B59-09ab at present, we cannot 
assure from our data only that the clump will not fragment later. We think, however, that 
further fragmentation is unlikely because fragme ntation tend s to proceed quite rapidly. For 
instance, some numerical studies, like those of IBossi (120091 ). show that oblate cores with 
magnetic fields can form binaries in timescales of less than 2-4 tff. Also, the models of 
Price fc Batd ( 120091 ) show that fragmentation in a core with initial mass of 50 M s un and 
initial radius of 0.375 pc -possibly not that much different from a low mass star cluster 
forming core like B59- proceeds within 1-2 tff independently of the amount of magnetic and 
radiative feedback support added to their models. 

Our analysis confirms the hypothesis that B59-09ab does not show significant fragmen- 
tation at the present time. We speculate that the clump is likely on its way to collapse, but 
it will not form multiple sources to increase the population of the small stellar cluster. The 
efficiency of formation in the cluster B59 after the collapse of the B59-09 clump will increase 
only modestly. B59 probably will remain as a small, low mas s star cluster with too few stars 
to survive disintegration by evaporation (ILada fc Ladall2003l ). 
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Fig. 1. — Top Left: MAMBO-II dust emission map of Barnard 59 from our observations. 
For purposes of clarity, the flux scale has been clipped from mJy/beam to 100 mJy/beam. 
Numbers refer to YSO sources identified and listed in Table [TJ The circle at the bottom 
right corner indicates the beam size. Top Right: MAMBO II Dust extinction map after 
conversion of the flux and subtraction of the YSOs emission. The extinction level has been 
also clipped from mag to 100 mag. The circle at the bottom right corner indicates the 
beam size. Bottom: NICER dust extinction map, as in RLA09, also scaled from mag to 
100 mag. The circle at the bottom right corner shows the size of the Gaussian filter used 
to construct the map. In all panels, the cross symbol marks the center around which radial 
profiles were constructed, at J2000 {a, 5) = (17 : 11 : 23,0, -27 : 25 : 59.3). 
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Fig. 2. — The open circle symbols show the radial profile of the NICER extinction map. 
The black diamond symbols show the radial profile of the MAMBO-II dust emission map. 
The profiles were constructed by averaging flux in circular, concentric rings centered on 
J2000 (a, 5) = (17 : 11 : 23,0, -27 : 25 : 59.3). For the MAMBO-II profile, the data were 
convolved to the resolution of the NICER map (20") and re-gridded to match pixels in both 
maps. 



Table 1. Properties of YSOs counterparts 



Gaussian fit 



lD a 


Class a 


Spectral Type b 


Peak Flux 


Mass b 


Aa 


A5 C 


deconvolved size d 








[mJy/beam] 


[M ] 


(" 


) 


("),(deg) 


BHB07-11 


1 




270.0 




-11.0, 


55.5 


18.0x17.0, -68.0 


BHB07-10 


0/1 




65.0 




65.0, 


-24.0 


16.0x16.0, 0.0 


BHB07-09 


Flat 


K5 


115.0 


0.77-0.79 


-29.4, 


-133.0 


13.7x12.1, -89.5 


BHB07-06 


II 


M2 


15.0 


0.24-0.62 


15.0, 


-80.0 


14.0x12.0, 79.4 


BHB07-07 


Flat 


K5 


15.0 


0.75-1.16 


15.0, 


-80.0 


14.0x12.0, 79.4 



a from iBrooke et all (120071 ) 
b from ICovev et all fcoioh 



c offsets from center of map at (a, <5) = (17:11:24.0, -27:25:30.0) 
d indicates major and minor axis, and position angle 
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Table 2. Barnard 59. Dust emission map. 



Parameter 


Value 


Tk a 


10.25 K 


rms 


4 mJy/beam 


total flux 


8.28 Jy 


peak flux 


90.7 mJy/beam 


Diameter b 


0.11 pc 


N ff2 c 


2.96 xlO 22 cm" 2 




1.30xl0 5 cm- 3 


Mass 


9.19 M 



Corresponds to our best fit 
to the Ay profile, not to a mea- 
sured value 

b Size of region with emission 
above 7„>0.2 I u , max 

c Average value over re- 
gion with emission above 



